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We investigate Big Bang nucleosynthesis (hereafter, BBN) in a cosmic environment characterised 
by a distribution of small-scale matter-antimatter domains. Production of antimatter domains in 
a baryo-asymmetric universe is predicted in some electroweak baryogenesis scenarios. We find that 
cosmic antimatter domains of size exceeding the neutron-diffusion length at temperature T ~ 1 MeV 
significantly affect the light-element production. Annihilation of antimatter preferentially occurs on 
neutrons such that antimatter domains may yield a reduction of the 4 He abundance relative to 
a standard BBN scenario. In the limiting case, all neutrons will be removed before the onset of 
light-element production, and a universe with net baryon number but without production of light 
elements results. In general, antimatter domains spoil agreement between BBN abundance yields 
and observationally inferred primordial abundances limits which allows us to derive limits on their 
presence in the early universe. However, if only small amounts of antimatter are present, BBN with 
low deuterium and low 4 He, as seemingly favored by current observational data, is possible. 
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Big Bang Nucleosynthesis is one of the furthest back- 
reaching cosmological probes available. By means of 
comparing the predicted and observationally inferred 
light element abundances the cosmic conditions as early 
as a few seconds after the Big Bang may be scrutinized. 
Nucleosynthesis has thus been used to constrain, for ex- 
ample, inhomogeneities in the baryon-to-photon ratio or 
abundances and properties of decaying particles during 
the BBN era (for reviews on non-standard BBN see jij). 

In this Letter we examine the BBN process in a uni- 
verse where matter and antimatter are segregated. A seg- 
regation of matter and antimatter may affect BBN abun- 
dance yields even if the average segregation scale is very 
small in terms of characteristic astrophysical scales. Such 
small scales imply annihilation of antimatter well before 
the present epoch. We are therefore particularly inter- 
ested in the case where the universe contains net baryon 
number, as for example, through an excess of matter 
domains over antimatter domains. A baryo-asymmetric 
universe filled with a distribution of small-scale matter- 
antimatter domains may arise during an epoch of baryo- 
genesis at the electroweak scale. It has been shown within 
the minimal supersymmetric standard model, and under 
the assumption of explicit as well as spontaneous CP vi- 
olation, that during a first-order electroweak phase tran- 
sition the baryogenesis process may result in individual 
bubbles containing either net baryon number, or net anti- 
baryon number Q. Recently, it has been argued that 
pre-existing stochastic (hyper)magnetic fields in the early 
universe, in conjunction with an era of electroweak baryo- 
genesis, may cause the production of regions containing 
either matter or antimatter J3j. In general, any segrega- 
tion of matter and antimatter either produced on scales 
larger than the neutron diffusion length at weak freeze- 
out, which corresponds to scales entering the horizon at 
cosmic temperatures below T « 1 TeV [0, or produced 



on superhorizon scales during earlier inflationary epochs 
may be subject to constraint by BBN. 

Nevertheless, the study of matter-antimatter domains 
present during the BBN era has to our knowledge not 
received any prior attention in this context. Earlier 
work on this problem Q has concentrated on the pres- 
ence of matter-antimatter domains in baryo-symmetric 
universes (i.e. Omnes-model). In essence, this work 
showed that antimatter domains and successful BBN mu- 
tually exclude each other. Another aspect, namely ho- 
mogeneous injection of antiprotons into the primordial 
plasma after the completion of a conventional BBN era 
(T < 10 keV) has been investigated by several authors 
g . The conclusion of such studies is that only a small 
fraction of about ~ a few 10~ 3 antiprotons per proton 
may be injected after the BBN epoch. Larger fractions 
lead to overproduction of mainly 3 He via p 4 He annihila- 
tion reactions. 

Annihilation of antimatter before, or during, the BBN 
era may be an astonishingly quiescent process. For lo- 
cal net antibaryon-to-photon ratios of the order, or less, 
than the net cosmic average baryon-to-photon ratio, 77, 
the entropy release associated with annihilation is only a 
minute fraction of the entropy in the photon-pair plasma. 
Heating of the plasma is therefore unimportant. In the 
absence of turbulent fluid motions and convective pro- 
cesses mixing of matter and antimatter should mainly 
proceed via (anti-) baryon diffusion. Neutrons diffuse 
most easily through the plasma and their diffusion length 
d n at the epoch of weak freeze-out (TV ~ 1 MeV) is ap- 
proximately d n « 10 4 cm. This corresponds to a comov- 
ing scale of 0.1 cm at T w 100 GeV; in what follows we 
often quote length scales comoving to this temperature 
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We distinguish three different limiting cases of the in- 
fluence of antimatter domains on BBN, according to the 
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segregation scale of matter and antimatter, or equiva- 
lently, the approximate matter-antimatter annihilation 
time: (a) Annihilations before weak freeze-out: No ef- 
fect on the light -element abundances because the n/p- 
ratio, which is important for the determination of light- 
element abundances, is reset to the equilibrium value by 
the fast weak interactions, (b) Annihilations between 
weak freeze-out (T « 1 MeV) and the onset of signif- 
icant 4 He synthesis (T rj 0.08 MeV): The n/p-ratio is 
strongly affected and a highly non-standard BBN sce- 
nario results. For reference, the neutron diffusion length 
scale at T a 0.08 MeV is roughly d n w 5 x 10 7 cm corre- 
sponding to a comoving scale of ~ 20 cm at T « 100 GeV. 
(c) Annihilations after 4 He synthesis (T < 80 keV): Pro- 
duction of elements and antielements in their respective 
domains occurs. Subsequent annihilation of antimatter 
proceeds mainly via p and p diffusion since n{n) are in- 
corporated into light elements. Due to the small mov- 
ability of p{p) in the plasma, this occurs mainly after the 
completion of BBN (T < 10 keV). The subsequent re- 
processing of light element abundances by annihilations 
may substantially deviate We are here most interested in 
the non-standard BBN scenario of case (b) . 

Annihilation of an antinucleon on a nucleon typically 
results in the production of a few pions and photons, 
i.e. NN i — ► afewTr^, 7T°, 7's. Pions may subsequently 
decay via 7r~ >—> fj,~ + ► e~ + v e + 9^ + and 
7T° 1— > 27. It is of interest if these pions or their decay 
products may alter the abundance yields, either through 
their effect on weak freeze-out or by, for example, the 
photodisintegration of nuclei. We assume that the num- 
ber of annihilations per photon is not much larger than 
a few x 10 -10 . In this case, annihilation generated v e 's 
have negligible effect on the 77,/p-ratio since their number 
density is more than nine orders of magnitude smaller 
than the thermal v e 's, which govern the weak equilib- 
rium. The same holds for electrons and positrons pro- 
duced in /i-decay which are quickly thermalized by their 
electromagnetic interactions. Annihilation-generated 7- 
rays cascade on the background photons via pair pro- 
duction and inverse Compton scattering, on a time scale 
rapid compared to the time scale for photodisintegration 
of nuclei [|o|. This cascade only terminates when in- 
dividual photons do not have enough energy to further 
pair-produce on background photons. For temperatures 
T > 3 keV, the energy of 7-rays below the threshold for 
e^-production does not suffice for the photodisintegra- 
tion of nuclei. Pions may perturb the n/p-ratio via charge 
exchange reactions, tt + +n t— > ir°+p and ir~ +p Tr°+n, 
provided they do not decay before. A typical inverse sur- 
vival time for pions against charge exchange reactions, 
(Un^dn^/dt « 5 x 10 5 (775/3 x 10~ 10 ) (T/MeV) 3 s -1 
jyj, should be compared to the decay rate of charged 
pions r^i w 3.8 x 10 7 s~ 4 , yielding a fraction / ~ 
0.01 (r) b /3 x lO" 10 ) (T/MeV) 3 of pions which may induce 
n 1— > p interconversion. Here r\\, is the local baryon-to- 



photon ratio in the annihilation region and it is under- 
stood that / is not to exceed unity. Assuming the produc- 
tion of about one pion per annihilation and 775 m 3 x 10~ 10 
it is straightforward to show that charge exchange reac- 
tions may affect Y p by no more than AY p < 2fR. This 
is small for small R and we will henceforth neglect the 
influence of pions on BBN abundance yields. 

The main effect of the presence of antimatter on BBN 
yields for annihilation between 1 MeV > T > 0.08 MeV 
and small antimatter fractions, R, is through a prefer- 
ential annihilation of antimatter on neutrons. Note that 
this effect had already been noted in the seventies || in 
the context of Omnes' baryosymmetric cosmology with 
the conclusion that such scenarios may conflict with a 
successful light-element synthesis. Neutrons diffuse out 
into the antimatter regions and annihilate there, whereas 
protons are confined to the matter regions. This results 
in a n/p-ratio reduction. Antineutrons which diffuse into 
the matter regions and annihilate on neutrons, as well 
as on protons, have a comparatively weaker effect on the 
71/p-ratio, if one assumes the density of antineutrons to 
be smaller than the density of neutrons. If annihilation 
occurs below a temperature of T < 1 MeV, weak in- 
teractions are nearly frozen out and the residual interac- 
tions cannot compensate for this annihilation-generated 
n/p-ratio reduction. Since the final 4 He mass fraction 
Y p strongly depends on the n/p-ratio at T rj 80 keV, 
in particular, Y p rj 2(n/p)/(n/p + 1), less 4 He than in 
standard BBN is produced in scenarios with antimat- 
ter. Abundance yields of other light isotopes are affected 
comparatively weaker. 

We have performed numerical simulations of the light- 
element nucleosynthesis in a universe with matter and 
antimatter domains. We have incorporated antinucle- 
ons and their weak interactions into the inhomogeneous 
nucleosynthesis code of This code solves the BBN 
network coupled to all relevant hydrodynamic processes, 
such as diffusion of baryons, photon diffusive heat trans- 
port, neutrino heat transport and late-time hydrody- 
namic expansion of high-density regions. The nuclear 
reaction network and the thermodynamics are treated as 
in the Kawano-code Jl3| . Note that antielements are not 
included in our code. 

Initial conditions for the determination of abundance 
yields are chosen to approximate a distribution of small, 
spherical antimatter domains immersed within ordinary 
matter regions. This is accomplished by performing sim- 
ulations with radial symmetry of antimatter bubbles sur- 
rounded by a comparatively larger matter shell. We have 
found that due to the large annihilation rates most anni- 
hilation takes place at the boundary between matter and 
antimatter domains Jlif . In our model, we have assumed 
constant and equal densities in the matter shells and an- 
timatter bubbles, r/i, — 7/5. Since we used the same cosmic 
spatial average matter density (77) = (1 — fv)f]b ~ fvVb 
in all simulations, this choice defines a unique relation 
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between the fraction of space in antimatter bubbles, fv, 
and the cosmic average ratio of antinucleons to nucleons 
R = fvVb/ {^—fv)Vb = fv/(l—fv)- It is understood that 
R is taken well after the QCD-transition so that there is 
no thermally produced baryon number, in particular, the 
77's reflect net (anti-) baryon-to-photon ratios. 

In Figure [j] we show the abundance yields of such 
matter-antimatter BBN scenarios as a function of the 
radius of the antimatter domains, ta, and for various an- 
timatter fractions, R. All calculations are performed for 
average net (rj) = 3.43 x 1CP 10 , corresponding to a frac- 
tional contribution of baryons to the critical density of 
il b = 0.0125/i -2 , with h the Hubble constant in units of 
100 km s" 1 Mpc" 1 . 
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FIG. 1. Mass fraction of 4 He and number densities of deu- 
terium and 7 Li relative to hydrogen as a function of the radius 
of the antimatter region r a in comoving units Q . Results are 
shown for different antimatter fractions R. 

It is evident from the top panel of Figure |l| that the 
presence of antimatter during BBN results in a reduction 
of the He abundance compared to the abundance in a 
standard BBN scenario (the values at the far left of Fig- 
ure 0) . The effect of 4 He-reduction is more pronounced 
for larger matter/ antimatter separation length scales. In- 
terconversion of neutrons into protons through residual 
weak interactions after the approximate weak freeze-out 
(Tyy) and neutron decay still reduce somewhat the n/p- 
ratio between weak freeze-out and synthesis of 4 He. The 
rates for these processes arc proportional to the neutron 
density. One may approximately assume that a fixed 
number of neutrons are removed through annihilation on 
antimatter for models with the same R. Early annihi- 
lation then weakens the post-weak freeze-out reduction 
of neutron density compared to annihilation taking place 
late. When segregation length scales become too large 
some annihilation of antimatter occurs after the bind- 
ing of neutrons into He. Since this is not yet properly 
treated in our code, the results on length scales larger 
than w a few cm in Figure ^ may be subject to modifi- 



cation, cf. pj. 

Within observational primordial abundance determi- 
nation uncertainties the abundances of the other light 
isotopes are only significantly modified with respect to 
standard BBN for large antimatter fractions. However, 
such large antimatter fractions produce observationally 
disallowed low 4 He mass fractions. It is interesting to 
note that there is no limit to the amount of 4 He reduc- 
tion possible. In particular, we have found scenarios with 
large antimatter fractions where all neutrons will be re- 
moved at temperatures above T » 80 keV. In such sce- 
narios there is also negligible production of other light el- 
ements. These would result in universes, or local regions 
of the universe, where no light-element production in the 
early universe occurs, but which, nevertheless, have net 
baryon number (cf. ||). 

In Figure || we compare approximate limits on the 
amount of antimatter from our BBN studies to similar 
limits derived from distortions of the cosmic microwave 
background radiation (CMBR). Whereas for ta between 
the neutron diffusion length at weak freeze-out and on- 
set of He synthesis the constraints derive from an un- 
derabundant production of 4 He, for larger ta limits are 
due to CMBR distortions. Input of non-thermal energy 
into the CMBR will be fully thermalized by means of 
Compton and double Compton scattering at redshifts 
larger than z — 10 7 . For the computation of the CMBR 
limits a characteristic annihilation redshift z of anti- 
matter domains of size ta was computed by the equal- 
ity, ta = d p (z), where d p (z) is the proton diffusion 
length. Limits were then determined by requiring that 
the amount of energy released by annihilations into the 
CMBR may not exceed the redshift-dependent observa- 
tionally allowed release of nonthermal energy into the 
CMBR as given in 

The BBN constraints on antimatter may be used, in 
turn, to constrain processes occurring in the early uni- 
verse before the BBN epoch. In the electroweak baryoge- 
nesis scenario by Q individual bubbles of broken phase 
during a first-order electroweak transition create cither 
net baryon or antibaryon number. Segregation of matter 
and antimatter is expected on the electroweak bubbles 
scale, Ijf, which has been estimated to fall in the range 
1§ W ~ (10" 5 - lO- 1 )^, with 1% W the electroweak 
horizon scale 1§ W w 2cm(T/100GeV)- 2 (if). Observ- 
able signatures in the BBN yields for such segregation 
scales are not anticipated. Nevertheless, some scenar- 
ios of j!| envision only a slight excess in the number of 
bubbles, Nb, over the number of antibubbles, N b) i.e. 
(N b - N b )/(N b + N b ) ~ 10~ 3 . In this case the effec- 
tive segregation scale of a realistic distribution may well 
exceed the bubble scale. Meaningful limits on the mod- 
els by H from BBN are more probable. Here stochas- 
tic matter-antimatter distributions are produced by pre- 
existing (hyper)magnetic fields. Segregation scales are 
not related to l^ w , but to the coherence scale of the mag- 
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netic fields. Assuming equipartition between the energy 
in the magnetic fields and radiation, a scale-invariant 
magnetic field distribution, for example, may yield aver- 
age (anti-) baryon-to-photon ratios as large as |ry| ~ 10~ 5 
on the scale d n (Tw). Even though scale- invariance may 
be unrealistic, a similar scenario with a factor 10 5 smaller 
fluctuation amplitude already leaves an imprint on the 
BBN abundances. In any case, fluctuation amplitudes as 
large as |^| ~ 10 -5 have not been properly investigated 
in this Letter due to the influence of pions. It is further 
difficult to assess and model a realistic three-dimensional 
matter-antimatter distribution and its BBN yields. 
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FIG. 2. Schematic limits from BBN and CMBR distortions 
on the ratio of antinucleons to nucleons, R, allowed to be 
present in the early universe, rA is again given in comov- 
ing units Q. The constraint for large antimatter fractions is 
somewhat uncertain due to the neglected influence of pions, 
as discussed in the text. 

It is tempting to speculate about BBN yields from 
stochastic distributions with net cosmic r/ larger than 
that inferred from standard BBN. Local regions with 
large R do not contribute to the light element nucle- 
osynthesis but contribute to fib- If the protons of such 
regions are mixed with the yields of regions at low R 
and low as well as high 77, some relaxation of the bound 
on Qb inferred from standard BBN may result. It seems, 
however, not easy to satisfy simultaneously observational 
constraints on all the light-clement abundances at much 
larger fib. It is also not clear if such scenarios would 
require extreme fine-tuning of fluctuation parameters. 

In summary, we have investigated BBN in the presence 
of matter-antimatter domains. We have concentrated on 
the study of BBN with small antimatter fractions at the 
canonical cosmic net baryon-to-photon ratio ~ 3 x 10 -10 . 
Surprisingly, antimatter fractions of the order R < 1% 
surviving to the epoch of weak freeze-out may consid- 
erably alleviate the long-standing tension between the 
primordial He and deuterium abundances in standard 
BBN and their observationally inferred limits [Il8| . Even 



though other modifications of BBN have been proposed 
which also yield low 4 He and deuterium p9| , our model 
may have some motivation from electroweak baryogenesis 
scenarios which may generate antimatter domains. 
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